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Abstract 

We study the emission of fermion and gauge boson degrees of freedom on the brane by 
a rotating higher-dimensional black hole. Using matching techniques, for the near-horizon 
and far-field regime solutions, we solve analytically the corresponding field equations 
of motion. From this, we derive analytical results for the absorption probabilities and 
Hawking radiation emission rates, in the low-energy and low-rotation case, for both 
species of fields. We produce plots of these, comparing them to existing exact numerical 
results with very good agreement. We also study the total absorption cross-section and 
demonstrate that, as in the non-rotating case, it has a different behaviour for fermions 
and gauge bosons in the low-energy limit, while it follows a universal behaviour - reaching 
a constant, spin- independent, asymptotic value - in the high-energy regime. 



1 Introduction 



The leading motivation for studying higher-dimensional theories [HE] is that they provide 
a framework for the unification of gravitation with the rest of the fundamental forces. 
Gravity, and possibly scalar fields, propagate in a (4 + n)-dimensional spacetime (the 
Bulk), while ordinary matter is confined in a four-dimensional hypersurface, the Brane. 
In models with large extra dimensions, the traditional Planck scale is an effective scale, 
while the fundamental Planck scale of gravitation is related to it in terms of the number 
and size of the extra dimensions. In such a framework higher-dimensional black holes 
can be created in trans-planckian collisions [3]. This opens the possibility of seeing them 
in ground-based colliders [1] , or in cosmic ray interactions [5] . These higher-dimensional 
black holes and their properties have been the subject of a number of articles in the last 
few years - for a summary of their properties and phenomenological implications, see 

A black hole created in a high energy collision is expected to evaporate through 
Hawking radiation p] both in the bulk, through the emission of gravitons and scalar 
fields, and on the brane through the emission of fermions and gauge bosons. The black 
hole is expected to undergo a number of phases: First, the balding phase, in which 
the black hole emits mainly gravitational radiation, losing all "hair" inherited from the 
original particles. Then, in the spin- down phase, the black hole loses all its angular mo- 
mentum, through the emission of Hawking radiation. Third is the Schwarzschild phase, 
in which the black hole loses its actual mass due to the emission of Hawking radia- 
tion. Finally, the Planck phase - where the black hole's mass or temperature reach the 
characteristic scale of gravity - which needs a quantum gravity theory in order to be 
studied. There have been both numerical and analytical studies of the Hawking radi- 
ation emitted from such a higher- dimensional black hole, for the Schwarzschild phase 
[ini IIIl [I2l [la 112 [13 [Tg ini [III [13 120^ as wen as for the spin-down phase 

[2a[2l[23[26l[27l[28l[29l[30l[ni[3a[3^ 

In the present work we focus on the evaporation of a rotating (4 + n)-dimensional 
black hole through Hawking radiation in the form of fermions and gauge bosons on the 
brane. Numerical studies of fermion and gauge boson emission by a rotating black hole 
on the brane already exist in the literature [2S1 [221 [30] • In this article, we continue our 
effort, initiated in [31] with the study of scalar fields, to derive complementary analytic 
results for higher-spin fields using an analytical approach. In section 2, we consider the 
gravitational background corresponding to a higher- dimensional rotating black hole, and 
write down the equations, describing all radiative components of fermion and gauge 
boson fields propagating on the brane, in the form of a single master equation. In section 
3, we solve this equation analytically by using a well-known solution matching technique: 
we first derive the solution at the near horizon regime, and then at the far field regime. 
Next, we stretch and match the two solutions at the intermediate zone, in the low-energy 
and low-rotation limit, thus producing a smooth solution for all spacetime. We then use 
this solution in order to compute, for every type of field, the absorption probability, a 
quantity that characterizes the Hawking radiation. In section 4 we produce plots for the 
absorption coefficient and compare them with existing numerical results. In sections 5 
and 6, we study the asymptotic behaviour of the corresponding cross-section and the 
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profile of tlie energy emission rates, respectively. In section 7, we state our conclusions. 



2 Master equation on a brane embedded in a rotat- 
ing n)D black hole background 

The background around a (4 + n) -dimensional rotating black hole is given by the Myers- 
Perry solution [35j. As is usual, here we focus on the case where the black hole, being 
created by the collision of brane-localised particles, is characterised by only one non-zero 
angular momentum component parallel to our brane. Then, the line-element takes the 
form 

ds'=(l- -J^] de + ^^/i^ dtd^-^ dr' - E d0' 



where 



_ |^^2 ^ ^2 ^ ""'^^t-? ) ^ - ^Q^' ^ (1) 

A = r2 + a2--^, S = r2 + a2 cos^^, (2) 

and dQ"^ is the line-element on a unit n-sphere. The mass Mbh and angular momentum 
J of the black hole are then given in terms of the parameters a, /i 

Mbh = T-^-T^ , J = — Mbh a , (3) 

IottG n + 2 

with G being the (4-|-?T,)-dimensional Newton's constant, and An+2 the area of a (n + 2)- 
dimensional unit sphere given by 

2^{n+3)/2 

= r[(n + 3)/2] • 

The corresponding line-element on the brane can be found by projecting out the angular 
variables that parametrize the extra dimensions. In that case, the factor (if2^ disappears 
and the 4-dimensional brane background is described by the line-element 

fi \ . , 2afism^e ^^^^^ E . ^ . 



ds'=[l- 1 dt' + dtd^ - A - ^ 



r^ + a^ + ^^|^)sin^^V. (5) 



In the above < ip < 2tt, < 6 < tt and n stands for the number of extra, spacelike 
dimensions that exist transverse to the brane [D = 4 + n). 

We would like to study the emission of gauge bosons and fermions by the afore- 
mentioned projected black-hole background. We assume that the emitted particle modes 
couple only minimally to the gravitational background and have no other interactions. 
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therefore, they satisfy the corresponding free equations of motion. By using the Newman- 
Penrose formahsm [511 EZ] ? and assuming the factorized ansatz 



<fsit, r, e, if) = e--' e^"^^ i?.(r) S::^ie) , (6) 

where 5"™ (6*) are the so-called spin- weighted spheroidal harmonics [38l [39| BO f HT | H2l B3] . 
the free equations of motion for particles with spin s = 0, ^ and 1 may be combined to 
form the following complete "master" equation [29], satisfied by the radial part of all 
radiative components of the field. 



A-^ 4- ( A^+^ ^ 1 + 
dr \ dr 



7^2 _ igXA' 

+ iisur + s (A" - 2) 6s,\s\ - Asj 



Rs = 0, (7) 



where 

K = (r^ + a^) — am , Agj = Xgj + a^tu'^ — 2amuj . (8) 

In the above, Xgj is the angular eigenvalue appearing in the equation satisfied by the 
spheroidal harmonics, namely 

d ( „dS^.:\ / 9 9 9 „ „ (m + s COS 9)'^ , , 

sin^— ^ + aVcos2 0-2ac<;scos^- ^ . ' + A„- + g S?, = 0. 



sin 9 de \ d9 J \ sin^ 9 

(9) 

From Eq. ([7]), it is clear that the radial parts of the radiative components with s = |s| 
and s = — \s\ satisfy different equations due to the presence of the term that is multiplied 
by the 6s^\s\ factor. This obstacle can be overcome by redefining 

= A-I^lp+I,|, = P_|,| . (10) 

In terms of the new radial functions, the master radial equation on the brane takes the 
simplified form 



AN ( Ai-N ^]^{ ^ + 4..c.r - A, P.(r) = , (11) 



dr\ dr J \ A 

where the A"-term has now disappeared, and Asj = A\s\j + 2|s|. The angular eigenvalue 
A|s|j cannot be expressed in closed form, however it can be written in terms of a power 
series with respect to auj [381 HOI HH SS] as 



\s\j = -\s\i\s\ + 1) + V /fc (aiu)'' = jU + 1) - \s\i\s\ + 1) - -TT—rao; + ... . (12) 



By solving Eq. ( ITTl) one can compute the absorption coefficient for the prop- 

agation of the field on the specific gravitational background. This quantity is needed 
to compute the Hawking radiation of the black hole on the brane. For example, the 
differential energy emission rate is given by 



dt duj 27r exp Ik/Tij] ± 1 
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with k and Tjj given by 

,_ o "'^^ nr {n + 1) + {n - l)al 

k = u-mQ = u--^^, Tn = ^^^^^-^^. (14) 

with VL the angular velocity, Th the temperature of the black hole and th the black-hole 
horizon radius, defined through the relation A(rj:/) = 0. 



3 Analytical Solution 

In this section, we will derive analytic solutions to the radial master equation valid at 
the two asymptotic regimes of the black- hole horizon (r ^ th) and far-field (r ^ th)- 
We will then demand that the two solutions are smoothly connected at an intermediate 
radial zone in order to construct a complete solution, valid at all radial regimes, for a 
field with arbitrary spin s. 



3.1 The Near-Horizon Regime 



In order to bring Eq. ffTTl) in the form of a known differential equation, we apply the 
following transformation [34] 



r f{r) 



A(r) + — yu/r 



n-l 



+ a? 



J.2 _|_ g2 



dr -|- 



(15) 



where A{r) = (n + 1) + (n — 1) a'^/r'^. Then, we may rewrite Eq. (fTTl) near the horizon 
(r ~ vh) - keeping s in the master equation as an arbitrary parameter - as 



d'^P dP 



df 



df 



L Alf(l-f) Alil-f) 
where we have defined oj^ = ujth, a* = a/rn. Also A'^ = ^'{th) = A^, is now 

2\s\+n{l + al) Aal 



P,(r) = 



5, = 1 - |s| + 
while A* and K^, are given by 



A^ 



Al 



A^ = n + l + {n- l)al , 
K^ = {1 + al)uj^ — a^m . 



(16) 
(17) 

(18) 
(19) 



If we make the redefinition: P{f) = /"(I — f)^F{f), Eq. (fT6l) takes the form of a 
hypergeometric equation ^6] 



(20) 



d'^F dF 
f{l-f)^ + [c-{l + a + b) f] — -abF = 0. 
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with 

a = a + (3 + B^-l, b = a + (3 , c=l-|s|+2a. (21) 



The power coefficients a and f3 can be determined by solving the modified second-order 
algebraic equations 

a^-|.|a + %-^ = (22) 



and 



respectively. The first of these two equations leads to the following solutions for the 
parameter a 

\S\ , / iK^: S~ 



while the second equation for (3 admits the solutions 
1 



(25) 

The general solution of the master equation near the horizon is then given by 

PNH{f) = A^r{l-ffF{a,h,c-f) 

+ A+r"(l-/)'3F(a-c+l,6-c+l,2-c;/). (26) 

We must now impose the boundary condition that no outgoing modes exist near the 
horizon of the black hole. Using the solutions (IMl) . in the limit r — * th, or /(r) 0, we 
obtain either 

for a = a+, or 

PMf) ^ ^-/^ r""'^^' + A^r'-^ r""*'^* , (28) 

for a = a_. We may now use the tortoise-like coordinate y = r^il + al) ln(/)/74^, that, 
in the limit r ru, becomes identical [31] to the usual tortoise coordinate r^,, defined 
by dr^/dr = (r^ + a^)/A(r). Then, the factors f^'^^*/^* reduce to e^*'^^ describing an 
outgoing and incoming free wave, respectively. According to Teukolsky's classic analysis 
[17], the correct boundary condition at the horizon of the black hole, for a field with 
non-zero spin, is 

Rs ~ A-^ e"^^"* , (29) 

which, in our case, translates to 

P+\s\ ~ e-^'^^ = , P_|,| ~ f\'\ e-'^y = /'"I . (30) 

Demanding that our asymptotic near-horizon (NH) solution obeys the above boundary 
condition, leads to the selection a = a_ and Aj^ = 0, bringing the NH solution to the 
final form 

PNHif) = A_ r (1 - fy Fia, b, c; f) . (31) 
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The criterion for the convergence of the hypergeometric function -F(a, b, c; /), i.e. Re{c — 
a — b) > 0, can finally be applied, and leads to the choice /3 = All the above results 
reduce to the ones for scalar fields [M] if we set s = 0. They also reduce to the ones valid 
for general spin-s fields in a non-rotating black hole background [7] if we set a = 0. 

For the purpose of matching the near-horizon and far-field (FF) solutions at an 
intermediate radial zone, we need to extrapolate ('stretch') our NH solution to values of 
the radial coordinate that are much larger than the horizon radius. We will do that by 
changing first the argument of the hypergeometric function from / to 1 — / by using the 
following relation [16] 



PnhU) 



A-r (1 - fY 



T{c)T{c-a-b) 
T{c-a)T{c-b) 



F{a, b,a + b - c+ 1;1 - f) 



+ (1-/) 



c—a—b 



r(c)r(a + b-c) 



r(a)r(6) 

The function /(r) may be alternatively written as 

/i 1 



F(c — a, c — b, c — a — b + 1] 1 — f) 



.(32) 



1 - 



1 



r 



n-l 



(33) 



where we have used the horizon equation A(r/^) = in order to eliminate /i from the 
above relation. In the limit r ^ Th, and for n > 0, the above expression goes to unity. 

By using the above, the argument of the "stretched" hypergeometric function goes 
again to zero, and the "stretched" near-horizon solution takes the form 

P (f\r^A ^^ ^^/^ ^(c)^(c-a-&) ^^_^+2-b.-|.| T{c)T{a + b-c) 

PMf) ^ A_ {i-ff _ _ + A_ (1 -/) r(a)r(&) • (^^^ 

For r 3> rjf, the quantity (1 — /) can be accurately approximated by 



1-/ 



n+l 



thus, bringing Eq. ( !34l) to a simpler power-law form 

PNH{r) ~ Air-("+i)^ + A2^(«+i)(/3+kl+B.-2) 

with 



A, = A, [{1 + al) rl^'] 



13 r(c)r(c — a — b) 



V{c-a)V{c-b) ' 
(p+\s\+B,-2) r(c)r(a + b-c) 



T{a)T{b) 



(35) 

(36) 

(37) 
(38) 
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3.2 The Far-Field Regime 

We now turn to the far-field regime, where the radial master equation (1111) becomes 

d'P, ^ 2(1 -\s\) dP, ^f2^^_ Xw + 2\s\+a'u' \ ^ ^ ^^^^ 



dr'^ r dr \ r 

By redefining Pg = e~*'^^r5(2MI-i+-^)p^^ with 



Z = ^{2\s\ - 1)2 + 4(A|,|j- + 2|s| + a2cj2) , (40) 

Eq. (!39|) takes the form of a confluent hypergeometric differential equation whose solution 
can be expressed in terms of the Kummer functions M and U, i.e. 



Pppir) = e-^'^V5(2kl-i+^)x 

[5i M (I - s + f , 1 + Z, 2tujr) + ^2 t/ (I - s + f , 1 + Z, 2tujr) ] . (41) 

As in the case of the NH solution, the FF one also needs to be "stretched", this time 
towards small values of the radial coordinate. In order to do this, we take the r — 
limit in (HTj) . which gives [16] 



PMr) - B, ri(2kl-i+^) + B2 r^(2N-i-^)^Z(^) (2,^)-^ . (42) 

r(2 — 5+2") 

In analogy with the stretched NH solution fl36p . the stretched FF solution also has 
a power-law form. In order to construct a complete radial solution, we have to match 
these two expressions at an intermediate zone. To this end, we note that, for urn <^ 1 
and a* -C 1, we get, from Eq. f|T7|) . 5* ~ 2 — |s| + , and 



P ^ ^(^^ (l - 2k| - 7(2|.|-1)2 + 4A.,) . (43) 
In this case, the stretched near-horizon solution fl5B]) takes the form 

P^^(r) = ^^^K2|^|-W(2|«hl)^+4A.,) ^ ^^^i(2H-l-V(2N-l)H4A.,) ^^^^ 

From the definitions (IE]) and ( l40l) for the constants A\s\j and Z, respectively, one can see 
that, within our approximation, the powers of r in the two stretched solutions ( H2l) and 
(1441) become identical. Then, a smooth matching is realized provided that we identify 
the corresponding coefficients in front of the same power of r. From this, we obtain 

B2 r(i-s + f) A2 

with A1/A2 following from Eq. fl38|) 



Al _ r{c - a - b)r{a)r{b) r 2n n+11 2l3+\s\+B.-2 . ^ 

A2 ~ r(c - a)r(c - b)ria + 6 - c) + J • ^^^^ 

Equation fHS]) ensures the smooth matching of the two asymptotic solutions, and thus the 
existence of a complete solution to the radial master equation describing the propagation 
of a field with arbitrary spin s in the gravitational background induced on the brane. 
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3.3 Computing the Absorption Probability 



The derivation of the complete solution for the radial function Ps{r) opens the way for 
the calculation of the corresponding absorption probability lA^j^p. The latter appears 
in the expression of the various emission rates for the Hawking radiation emitted on the 
brane from the higher-dimensional black hole given in Eq. ([1]). In order to derive this 
quantity, we expand the far- field solution ( |4T1) in the limit r ^ oo, and obtain [46] 



^ , , , 5ie*n2--+2j T{l + Z) 

(2za;)(^^+#)r(i + 3 + f) (2. 




+ 



Bi r (1 + Z) e^'^^ 

_ ^ + (47) 

Let us focus first on the case of fields with spin 1/2: from the master radial equation 
TTl) one may easily derive, similarly to the 4-dimensional case [371 SHI HHl EO], that 

^(\P,\^-\P_,\A =0. (48) 



dr V 2 

The conserved - for arbitrary values of r - quantity inside the parenthesis is proportional 
to the radial component of the particle current produced by the black hole. Then, the 
absorption probability is defined as the ratio of the flux of particles at the black hole 
horizon over the one at infinity 

in in 

where in the second part of the equation, we have used the conservation of the total flux. 
The flux of fermions at infinity may be found by integrating the radial component of the 
conserved current over a 2-dimensional sphere at asymptotic infinity. Applying this and 
using Eq. P7|) . we find 

|p(o"*)|2 _ |p{o^"*)|2 |yoMt|2 

l^ljml —1 |pM|2_|pM|2 ~ ~'~ |yM|2 ' ^^^^ 

2 ~ 2 2 

Using, finally, the explicit expressions for y}°^^\ as these are defined in Eq. (H7|) . we find 

^2 

that they are related through the equation 

y^r^ = , (51) 

which, therefore, leads to the following expression for the absorption probability for 
fermions 



Ui™ = 1 



Y(out) 2 



(in) 
I 

2 



(52) 
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For fields with spin s = 1, there is no conserved particle current. In order to compute 
the absorption probabihty, a technique, introduced in |f5T|, may be followed in which the 
radial master equation is transformed, through a radial function redefinition and the use 
of the tortoise coordinate r^,, to an alternative one with real, short-range potential. Then, 
the asymptotic solution at infinity for the gauge field is given, in terms of the new radial 
function, by the following expression [SH E2] 



(53) 



i.e. by the sum of outgoing and incoming plane waves with constant amplitudes, from 
which the expression of Aij^ may easily follow 



14,. |2 

I ^Ijm I 



1 - \A 



(in) 1 2 



(54) 



Although the exact analysis is quite cumbersome [STl [52] . it yields relations connecting 
the constant amplitude A^j^^ in Eq. ( I53i) with the constant coefficients appearing 
in our Eq. (H71) and leading finally to 



I ^Ijm I 



1 - 



16a;4 

jmu) 



(out) 



(in) 



(55) 



The constant Bjmuj is defined as the coefficient appearing in the differential equation 

A dIdIp+i = Bjmu. P-i , (56) 
where = dr + iK/ A, or, equivalently, through the relation 



yiout) 



B 



(out) 



(57) 



holding for the asymptotic solution P7|) when substituted in Eq. fl56l) . By using the 
explicit expressions of from Eq. (H7|) . we obtain 



Y^out) 



(out) 



{51 



which leads to Bj^uj = Xij + 2 + a u , and, thus, to the final expression for the absorption 
probability for brane-localised fields with spin 1 



Ijm I 



{Xij + 2 + a^Lo^y 



(out) 



(in) 



(59) 



4 Plotting our Analytic Results 

We can now use Eqs. and flS^ - our main analytic results - to produce plots for the 
absorption probabilities for spin-1/2 and spin-1 particles, respectively, propagating in 
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|Ai|2 |Ai|2 

2 2 




Figure 1: Absorption probability |^i/2p for brane spinor particles, for the modes j = 1/2 
and m = 1/2, —1/2, from left to right, for n = 6 and a* = 0.0, 0.5, 1.0, 1.5. 



|A,|2 lAip |AiP 




Figure 2: Absorption probability \Ai\'^ for brane boson particles, for the modes j = 1 
and 171 = 0, —1, 1, from left to right, for n = 6 and a* = 0.0, 0.5, 1.0, 1.5. 



the induced-on-the-brane gravitational background. Upon selection of particular values 
of the angular momentum numbers {j, m), the values of Agj^ will be plotted as a function 
of the energy parameter urn, and then compared with existing numerical results from 
the literature [281 [29l [30] . 

We remind the reader that in the process of producing a complete analytic solution 
to the radial master equation, we were forced to assume that uth <^ 1 and a* -C 1, 
therefore, strictly speaking, our results are valid only in the low-energy and low-rotation 
limit. However, at times, our plots will extend beyond this range of validity to exhibit the 
remarkably good agreement we obtain even outside the assumed range of validity of our 
approximations. Another point that we would like to stress is that, in general, an increase 
in the number of extra dimensions n improves the validity of our approximation: as n 
increases, the assumed behavior of /(r) at infinity in Eq. (JT5i) becomes more accurate, 
and terms that were neglected during the matching of the asymptotic solutions, such as 
K^:/A^, become even more suppressed. 

In Figs. [1] and [21 we plot the absorption probability for fermions and gauge bosons, 
respectively, for the lowest partial waves in each case, i.e. (j = 1/2, m = ±1/2) and 
(j = 1, m = 0, —1, 1), in the case of n = 6 and for various values of a^,. From Fig. [H we 
see that, for fermions with m > 0, the absorption probability is monotonically increasing 
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with the angular momentum parameter over most of the energy regime, while, for m < 0, 
|74i/2p increases with a^, in the low-energy regime and is suppressed at the high-energy 
regime. For bosons, and m = 0, the absorption probability is increasing with the angular 
momentum of the black hole for small values of a* but this increase is decelerated (and 
eventually reversed) for higher values of a*; for m < 0, the absorption probability is 
enhanced with a^, in the low-energy regime but suppressed in the high-energy one; finally 
for m > 0, an increase in a^, leads to more negative values of in the superradiant 
regime [S3] (where the quantity uo — rriQ, as well as l^ip, becomes negative signalling 
the enhancement of the incoming wave's amplitude by a rotating black hole) and to 
more positive values in the non-sup err adiant one. Our results are in excellent agreement 
with the exact numerical ones p8l [29l [30] in the low-energy and low-angular-momentum 
regime, and remain in very good agreement even beyond this range of values - a deviation 
from the exact results appears only for values of the angular momentum parameter larger 
than unity. 

A similar agreement with the exact behaviour appears in the dependence of the 
absorption probabilities for both fermions and gauge bosons on the number of extra 
dimensions n. To avoid repetition of known results, we do not present here any plots for 
the dependence of |y4i/2p and |y4ip on n, but only comment on the derived behaviour. For 
fermions, a monotonic decrease is observed for m < 0, while for m > 0, an enhancement 
in the low-energy regime is followed by a suppression in the high-energy one. For bosons, 
a monotonic suppression with n is found for the non-superradiant modes (m < 0) over 
the whole energy regime, while for the superradiant ones (m > 0) a similar suppression 
in the non-superradiant regime is replaced by an enhancement (i.e. more negative values) 
in the superradiant one. 

5 Asymptotic Values and Absorption Cross-Section 

By expanding further our analytic expressions (1521) and (!59|) for the absorption proba- 
bilities for fermions and gauge fields, respectively, one may obtain simplified, compact 
expressions that reveal more clearly the low-energy asymptotic behaviour in each case 
as well as potential differences in the asymptotic values of the corresponding absorption 
cross-sections. 

In the limit u 0, we obtain, from Eq. fl40l) . Z ~ 2j + 1, and then, from the 
expression of B1/B2, Eq. ( HSl) . the result 

f - (2...)--^ (^) ^ M.,^ {2^u;r^-' , (60) 

B2 r(i+j-s) v^2/^=o 

with Mgjm a complex constant independent of the energy u. Starting from the fermions 
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Figure 3: Total absorption cross-section for spinor fields, in rj^ units, for the cases (a) 
n = 6 and various a*, and (b) a* = 0.5 and various n, in low energies. 



and Eq. (152|) . we set s = 1/2 and find 



\A, I 



Mi,^r(2j + 2) 



Mi.T{2j + 2) e-0-+V2) + r(j + 3/2) (2icu) 

2 



2i+l 



(2cu)^^-+ir(j + 3/2) 

r(2j + 2) 



M 



1 1 

+ 



■jm 



Ml. 



+ 



The corresponding absorption cross-section, defined as agjm 
partial wave [5^, will then have the form 



(61) 



IT \Asjm\'^ /^^'^ for each 



_ 2^J'+%c^^^'-^r(j + 3/2) / 1 1 

- r(2j + 2) I mT" m 



+ ... . 



(62) 



From the above result, we may easily read that, similarly to the case of a non-rotating 
higher-dimensional black hole [10lll2j, the absorption cross-section of the lowest fermionic 
mode, with j = 1/2, assumes a non-zero asymptotic value, namely 



o-i 1, 

2 2 



27r 



Mil 



MX. 



(63) 



while all higher fermionic partial waves, with j > 1/2, will have zero absorption cross- 
section as — i> 0. The quantity Mii^ depends both on the number of extra dimensions n 
and on the angular momentum parameter a* of the black hole. In Figs. [3](a,b), we depict 
the dependence of cri - summed over j and m - on both a^, and n, respectively. We observe 
that the asymptotic value of the absorption cross-section for fermions in the low-energy 
regime is enhanced in terms of both parameters of the gravitational background. 

For gauge fields, with s = 1, Eq. (l59ll leads to a similar result for the absorption 
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Figure 4: Total absorption cross-section for vector fields, in rj^ units, for the cases (a) 
n = 6 and various a^,, and (b) = 0.5 and various n, in low energies. 



probability, namely 

Mi,„r(2j + 2) 



I -^Ijm 



2 



Mij^T{2j + 2) e-^- + r(j + 2) {2tuj)^^+^ 
(2^)2^-+ir(j + 2) z(Mr,„-Mi,J 



r(2j + 2) |Mi,^|2 + ^^^^ 

The corresponding absorption cross-section will similarly have the form 

In this case, all gauge fields modes, including the lowest one with j = 1, have zero 
asymptotic value of absorption cross-section - this is in agreement with the corresponding 
results derived in the non-rotating case [T0l[T2]. The dependence of cxi - summed again 
over j and m - on a* and n is now given in Figs. |l](a,b). Again, an increase in both 
parameters causes an enhancement in the absorption cross-section for gauge fields in the 
low-energy regime. 

As we saw in the previous section, the absorption probability l^sj^P for both species 
of fields, fermions and gauge bosons, approaches unity for large values of the energy 
parameter uurn- Therefore, the absorption cross-section for each partial wave will start 
from either a zero or non-zero low-energy value - depending on the spin of the field and 
its partial wave numbers, will reach a maximum value and then asymptote to zero as 
Oi^l/uj'^). Nevertheless, we expect the behaviour of the total cross-section cr^ - summed 
over j and m - to be radically different. Analyses in the case of a non- rotating black hole 
fl2\ [55] have demonstrated that the total cross-section for all brane-localised species 
of fields - scalar, spinor and gauge fields - reach asymptotically a universal constant 
high-energy value that depends only on the number of extra dimensions. In the case 
of a rotating black hole, it has been similarly shown, both in the 5-dimensional [33] 
and (4 + n)-dimensional case [31] that, in the case of scalar fields, a similar constant 
asymptotic value is reached in the high-energy regime. In the latter work [34], it was 
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Figure 5: Total cross section, in rfj units, for spinor and vector fields, for the case n = 6 
and a* = 0.5. The dashed line denotes the theoretical asymptotic value of cr, as computed 
in [31]. 



shown in detail that the absorption cross-section can be described in the high-energy 
regime by the geometrical optics limit and that the exact asymptotic value of ag can be 
approximated by an analytic expression giving the cross-section for a particle approaching 
the black hole parallel to the rotation axis. 

Thus, it remains to be seen whether the absorption cross-section, in the case of the 
remaining species of fields - i.e. fermions and gauge bosons, approaches again the same 
asymptotic value at high energies, as in the non-rotating case. By using our analytic 
expressions f l52l [59]) for the absorption probabilities, we have computed the partial wave 
cross-sections and, by taking the sum over j and m in each case, the total absorption 
cross-sections for fermions and gauge fields, respectively. Their behaviour over the whole 
energy regime is shown in Figs. [5](a,b). Although our results are expected to hold only 
in the low-energy limit and therefore not to yield any reliable results in the high-energy 
regime, we have nevertheless attempted to compute the high-energy asymptotic value 
of ag. Surprisingly enough, our results do not break down at high uvh and successfully 
yield a constant high-energy asymptotic value for the absorption cross-sections ai/2 and 
(Ti. In addition, the asymptotic values are identical to the one for scalar fields computed 
analytically in [34] - and denoted here with the horizontal dashed line. Our results, 
therefore, prove the universality in the behaviour of fields with arbitrary spin in the high- 
energy regime, and the dependence of their absorption cross-section only on parameters 
of the gravitational background even in the case of a higher-dimensional rotating black 
hole. 

6 Energy Emission Rates 

As is well known, the absorption probabilities lAsj^P can be used in order to determine 
the various emission rates for the Hawking radiation. In the present case, our derived 
analytic expressions ([S2D and (jSH]), together with the expression for the temperature ffn|) . 
can lead to the emission rates for a higher- dimensional rotating black hole directly on 
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Figure 6: Energy emission rates for spinor and vector fields, for = 0.5 and various n. 



our brane, in the form of fermions and gauge bosons. 




(orn 



Figure 7: Energy emission rates for spinor and vector fields, for n = Q and various a^,. 

Since our analytic expressions are, strictly speaking, valid in the low-energy and low- 
angular-momentum limit, we first focus our attention on these particular regimes. In the 
context of our analysis, we have studied all emission rates, associated with the Hawking 
radiation, namely the flux (number of particles), power (energy) and angular momentum 
emission rates, and they were all found to exhibit the same behaviour in terms of the 
number of extra dimensions n and the angular momentum parameter a^. Therefore, as 
an indicative example, we display here our results for the energy emission rates following 
from Eq. (fT3|) . In Figs. [6] and [71 we present the energy emission rate, i.e. the energy 
emitted by the black hole per unit time and unit frequency, in the form of fermions and 
gauge fields, as a function of n and a^, respectively. One may easily observe that the 
energy flux is indeed enhanced with both topological parameters, in accordance with the 
exact numerical results of [2Sl ED] • 

Even by looking at the low-energy regime, one can accurately conclude that the 
emission of gauge bosons dominates over the one for fermions, for the same values of 
the parameters a,, and n. This feature is also in accordance with the results of the exact 
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numerical analysis. In addition, the magnitudes of the corresponding energy emission 
rates are also correctly reproduced by our analytic results. Motivated by this, in Fig. [HI 
we plot the energy emission rates, for fermions and gauge bosons, for higher values of 
the energy parameter ujth and various configurations of ra, while a* is kept fixed. The 
restricted validity of our analytic results are bound to lead to deviations from the exact 
numerical results: indeed, our curves tend to reach their maximum point and their tail 
region sooner than the exact ones. Nevertheless, apart from the above, the agreement 
is remarkable: not only is the general profile of the energy emission rates and their 
dependence on the number of extra dimensions reproduced but the magnitudes at the 
peak of their curves - reached at ujth — 1, or even well beyond this - are also accurately 
derived. 



7 Conclusions 

In the present article, we studied the emission of Hawking radiation by higher- 
dimensional rotating black holes into fermionic and gauge boson degrees of freedom 
on the brane. The same topic has been studied before in the literature, however, those 
studies were either purely numerical or focused on specific cases like 5-dimensional space- 
times. In contrast, here, we carried out an analytical study of the problem deriving results 
valid for a black hole in a spacetime with arbitrary numbers of extra dimensions. 

As a starting point of our analysis, we formulated the radial master equation, de- 
scribing the propagation of degrees of freedom with spin s on the brane, in such a way 
that both radiative components of the given fields, i.e. with s = ±|s|, are described 
through the same equation. As the complexity of the problem forbids an all-energies an- 
alytical treatment, we followed an approximate method that allows us to find a solution 
in the low-energy and low-rotation regime. In Section 3, we used a well-known matching 
technique, which consists of solving the field equations in the near-horizon and far-field 
regime, and then matching the two solutions to construct a complete, smooth solution 
for the radial part of the field. This allowed us to compute the absorption probabilities 
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for brane-localised fermions and gauge bosons through two purely analytic expressions, 
namely Eqs. flS^ and flS^ . respectively. 

Using these expressions, the corresponding absorption probabilities were plotted as 
a function of the energy-parameter ujth, for various values of the number of extra di- 
mensions n and angular momentum parameter a*. These plots were presented in Section 
4. Comparing our results to numerical results existing in the literature, we found them 
to be in very good agreement. The dependence of |y4sjmP on both n and a* was cor- 
rectly reproduced as well as additional effects such as the superradiance in the case of 
gauge bosons. This agreement extended beyond the low-u; and low-a* regime, even for 
intermediate values of a* or values of ujth larger than unity. In section 5, we presented 
compact, simplified expressions for the absorption probabilities valid in the limit 0. 
As in the non-rotating case, the behaviour was found to be different with the absorp- 
tion probability for gauge bosons reducing to zero faster than the one for fermions. This 
resulted in a different behaviour, in the very low-energy regime, of the corresponding 
total absorption cross-sections: while the one for gauge bosons goes always to zero, the 
one for fermions assumes a non-zero, constant value that depends on both a* and n. 
The extended validity of our results allowed us to compute the dependence of the total 
cross-section over the whole energy regime. We were then able to derive the high-energy 
asymptotic values, that turned out to be the same for fermions and gauge bosons and to 
coincide with the one for scalar fields computed in a previous work of ours via the use 
of the geometrical optics limit. 

Finally, in section 6, we computed the energy emission rates for Hawking radiation by 
a rotating, higher-dimensional black hole on the brane in the form of fermions and gauge 
bosons. These were also found to be in remarkable agreement with the exact numerical 
results in the low-energy and low-angular-momentum regime, and to lead also to fairly 
accurate results beyond those limits. 

In conclusion, our study shows that analytical treatment of the propagation of fields 
with spin s on a brane embedded in a higher- dimensional, rotating black hole background 
is indeed possible, and that the derived results, although they cannot offer the level of 
accuracy and completeness of the exact numerical ones, can certainly reproduce quite 
accurately the dependence of the various quantities on the topological parameters of 
spacetime as well as a good estimate of their actual magnitudes. 
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